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The Darboux transformation operator technique is applied to construct 
exactly solvable anharmonic singular oscillator potentials and to study 
their coherent states. Classical system corresponding to a transformed 
quantum system is constructed with the help of the coherent states tech- 
nique. It is shown that at classical level the Darboux transformation may 
be treated as a transformation of Kahler potential which leads to a dis- 
tortion of the initial phase space. 
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I. INTRODUCTION 

In quantum mechanics there exists a beautiful recipe to construct a family of 
exactly solvable Hamiltonians starting from an initial Hamiltonian ho for which 
the solutions of the Schrodinger equation hoipE = EipE are known. Since this 
method is very simple in applications and possesses an intrinsic charm it attracts a 
thorough attention of several generations of researches. A central point of this theory 
is a notion of a transformation operator introduced and investigated by Delsart 
and Lions.0 According to their approach, every operator L intertwining the initial 
Hamiltonian ho and transformed Hamiltonian hi, 



is a transformation operator (originally called a transmutation operatorlH). It is 
obvious that L transforms solutions ipE {$e ^ kerL) of the initial Schrodinger equa- 
tion onto solutions tps — LipE of a transformed Schrodinger equation h\ifE = EtfE 

a To appear in J. Math. Phys. 1998, 39, No 1 or No 2. 
b Electronic mail address: samsonov@phys.tsu.tomsk.su 
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corresponding to the same eigenvalue E. In the case of a first order differential 
operator L and a one dimensional Hamiltonian ho, the relation (|I|) leads to! the 
transformation proposed by Moutardl and studied by Darboux @ (see also Ref.|5]). If 
L is a differential operator of order N it is called! N -order Darboux transformation 
operator. In this paper we will restrict ourselves by first order transformations. 

The Darboux transformations operator technique is a very useful tool in obtain- 
ing and investigating new exactly solvable quantum problems. It enters many fields 
of physics and mathematics. In soliton theory the Darboux transformation permits 
one to construct and investigate properties of multisoliton solutions to nonlinear 
evolution equations (see Ref.^| and references therein). In quantum mechanics this 
technique takes the form of a factorization method. Using this method Schrodinge r0 
was able to solve the hydrogen atom problem in a pure algebraic manner. Subse- 
quently, Infeld and Hali generalized the Schrodinger method and obtained a wide 
class of solvable potentials. In the inverse quantum scattering problem! this trans- 
formation is used to construct potentials with given properties and to obtain appro- 
priate solutions to the Schrodinger equation. New developments in this field have 



recently been presented in Ref.|10|. In the context of supersymmetric quantum me- 
chanics (for reviews see Ref.|TTD the Hamiltonian ho and its superpartner hi form 
a matrix-valued Hamiltonian (so called super Hamiltonian) TC. Darboux transfor- 
mation operator L and its Hermitian conjugate L + participate in construction of 
supercharge operators which, together with T~C, form a superalgebra. Initially the 
supersymmetric quantum mechanics served as a toy model to illustrate the problem 
of supersymmetry breaking in quantum field theories. Recently, new features in s 
upersymmetry breakdown in quantum mechanics caused by higher order Darboux 
transformations have been discovered.il Besides quantum field theories the same 
ideas find applications in many others fields of theoretical physics such as nuclear, 
a tomic, solid state and statistical physics.lD] 

In spite of the long period of application of the Darboux transformation operator 
method to quantum mechanical problems the question of correspondence between 
the transformed quantum system and a classical system has not been properly dis- 
cussed in the literature. In a recent paperll! a classical mechanics counterpart of 
the transformed system has been described. It consists in addition of a total time 
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derivative of a purely imaginary function to the Lagrangian of a classical system. 
This transformation does not affect the equations of motion. Canonical quantization 
of the new classical system gives the transformed quantum system. 

The concept of coherent states plays an important role in several fields of physics 
and mathematics (see, for example, Refs. p^ - |l5|) . In this paper we exploit the prop- 



erty of coherent states to establish the correspondence between the quantum me- 
chanics and the classical mechanics. For any quantum system we can construct 
using the coherent states technique!!' HJ a classical system possessing the following 
remarkable property. The geometric quantizationlll of the resulting system gives 
holomorphic representation of the initial quantum system. By these means we may 
establish a correspondence between quantum and classical systems. Recently it has 
been notedlll that application of the Darboux transformation operator to the coher- 
ent states of the initial quantum system gives a set of states which may be treated 
as coherent states of the transformed system. This observation open s the door to 
establishing a correspondence between transformed quantum and classical systems 
and obtaining a classical counterpart of the Darboux transformation. 

In this paper we realize the above program for a system with the Hamiltonian 

h = -d 2 /dx 2 + x 2 /A + b/x 2 , 6>0, a?e[0,oo). (2) 

We will not discuss the significance of the spectral problem for the Hamiltonian 
(0) (more precisely, for its closure in the Hilbert space H of square integrable on 
interval [0, oo) functions with the zero boundary condition at the origin). We only 
notice the fact that this Hamiltonian appears in different physical problems such as 
central potential problem, Calogero problemlll, fractional statistics and anyonsI3, 
two- dimensional qcdIH, and others. 

We obtain exactly solvable partners for the Hamiltonian (|j). Next we construct 
the coherent states for the transformed a Id Darboux Hamiltonians. These states as 
well as the coherent states for the input system (initial coherent states) are labeled 
by the points of a unite disk on the complex 2-plane. We consider two different 
Kahler structures on the unit disk. The first one is the well known Poincare model 
of the Lobachevsky plane and the second is the Darboux transformation of the first. 
The expression for the transformed Kahler potential is obtained. We also notice 
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that in contrast to the initial manifold the transformed one has a varying curvature. 
The integration measure realizing the identity decomposition for the transformed 
coherent states is calculated. The holomorphic representation for the discrete basis 
set as well as for the coherent states, lowering and raising operators and Darboux 
transformation operators is constructed. It is established that at the classical level 
the Darboux transformation may be considered as such a transformation of the 
Kahler potential and Hamilton function that preserves the equations of motion. 



II. COHERENT STATES OF SINGULAR OSCILLATOR 

In this section we briefly review the well knownEl properties of the coherent 
states for the Hamiltonian (0) that we will need further. 

The dynamical symmetry algebra for the quantum system with the Hamiltonian 
(0) is su(l,l). In coordinate representation the generators of this algebra are ex- 
pressed in terms of the harmonic oscillator annihilation a and creation a + operators 

k = i/i , k + = \ [(a+) 2 - b/x 2 ] , k_ = \[a 2 - b/x 2 } , 
a = d/dx + x/2, a + = — d/dx + x/2. 

and satisfy the standard commutation relations 

[h, k±] = ±k±, [k-, k+] = 2k . 

Since we consider an irreducible representation the corresponding su(l,l) Casimir 
operator takes the constant value C = ~ [k + k^ + k^k + ] — k 2 , = 3/16 — 6/4 = k(l — k). 
The value of /c = 1/2 + (1/4)a/1 + 46 defines the ground state (vacuum) energy 
E = 2k. The vacuum state |0) is defined by the equations: A;_ | 0) = and 
ko | 0) = k | 0). Other discrete basis eigenf unctions | n) of the Hamiltonian ho are 
defined by means of the raising operator k + 



k + | n ) = -yj(n + l)(n + 2k) \ n+ 1), 



n) = (-1)' 



li2/ "' ) (£:+)" | 0), k \n) = {k + n)\n) 



\\ n\T(2k + n 

and in the coordinate representation they have the form 

ip n (x) = (x | n) = [n!2 1 " 2fc r- 1 (n + 2k)] ^ x 2 ^ 2 exp(-x 2 /A)Lf-\x 2 /2), (3) 
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where L%(z) are the Laguerre polynomials. 

The coherent states | z) can be obtained by applying of the unitary displacement 
operator D z to the vacuum state 



z) = D z \0) = e zk+ exp [ln(l - zz)k ] e~ zk ~ | 0) = No, £ a n z n \ n) 

n=0 



N 0z = (1 - zz) k , a n = (-l) n t 



T(2k + n) 

\z\ < 1. 



\ n\T(2k) 

In the coordinate representation they look like 

fa(x) = 2 1 / 2 - fc r- 1 / 2 (2A;)(l - z)- 2k (l - zz) k x 2k -V 2 x 
exp [-(1/4)(1 - z)-\l + z)x 2 } , |*| < 1. 



(4) 



The linear manifold spanned by the vectors {| z)} forms an everywhere dense 
set in the Hilbert space H with the following identity decomposition: 

2k — 1 

I z)(z I dfj,(z) = 1, dfj,(z) = (1 — z~z)~ 2 dzdz. (5) 

|z|<l 7T 

The system {| z)} is overcomlete. The completeness of different subsystems of this 
system was studied in detail in Ref.|22|. 

The Fourier coefficients c n of a vector | tp) e H with respect to the basis {| n}} 
define a holomorphic representation ip(z) of the vector | ip) in the space of functions 
which are holomorphic in the unit disk 

oo 

(z\ip} = N 0z tp (z) , ijj (z) = 53 a n c n z n . 

n=0 

The functions ip n {z) = a n z n realize the holomorphic representation of the basis 
vectors | n). For the holomorphic representation of the coherent state vectors we 



have 



C (0) (z)=N - 1 (z\() = (l-(() k (l-(z) 



-2k 



Let £ be the lineal of entire analytic functions ip{z) defined in the unit disk such 
that 

i>(z) | 2 (1 - \z\ 2 ) 2k dfi(z) < oo. 



\z\<l 

If now we define an inner product in the lineal £ 



(Vi (z) | ifa(z)) = I e-/ ( °Vi (z) ^2 (z) dfi (z) = (fa | fa), (6) 



z <1 
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where (ipi | ^2) is the inner product in the Hilbert space H then the completion 
of £ with respect to this inner product gives the Hilbert space 7i. The function 
= In |(0 I z)\' 2 = ln(l — zz) 2k is the Kahler potential in the unit disk. The 
corresponding metric has the form 

ds = gdzdz, 9=^ = j—f. 

As usually the imaginary part of this metric defines a symplectic 2-form uo = 
—igdzAdz and consequently a Poisson bracket {Hi, H 2 } of the functions Hi and H 2 
defined on the unit disk. The unit disk represents now a phase space of a constant 
Gauss curvature IC^ = — f . Making use of the covariant Berezin symbolslll we 
establish a one-to-one correspondence between the classical observables and the 
holomorphic representation of the quantum mechanical operators. For example, the 
following classical observables: 

1 + zz 



Kq — {z \ko \ z) — k 



I — zz' 



K_ = (z\k-\z) = = , K + = (z \k+\ z) 



1 — i r \ — r / -i — 

1 — ZZ 1 — ZZ 

correspond to the generators of su(l, 1) algebra of quantum operators. They form 
the basis of the algebra isomorphic to su(l, 1) with respect to the introduced Poisson 
bracket. 

If we choose Kq as a classical Hamilton function then the evolution of the system 
is described by the Hamilton equations 

Z = {Z, Kq] = -iz, Z = {Z, K } = %Z, 

where the dot stands for the time derivative. 

Geometric quantization^ of the obtained classical system^ leads to the holomor- 
phic representation of the quantum system described above. The value h = {2k)^ 1 
plays the role of the Plank constant. 

To conclude this section we present the explicit expressions for the Bergman 
kernel which takes off the integration in the £-space 



5® (z,z') = Y,^n (z)i> n (z') = (l-ZZ>^ 21 
n=0 



c In our case this is a quantization on the Lobachevsky planelil which really is a Berezin! 
quantization on a Kahler manifold 
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and for the su(l, 1) generators 

k (z) = z^- + k, k + (z) = z 2 ^- + 2kz, k_ (z) = 

dz dz dz 

in holomorphic representation. 

III. DARBOUX TRANSFORMATION OF COHERENT STATES 

In a simplest caseH i the Darboux transformation operator is the first order 
differential operator of the form 

L = —L (x) + d/dx = —u\x)/u(x) + d/dx (7) 

where the prime denotes the derivative with respect to x. When acting on the 
solutions ip n (x) of the initial Schrodinger equation 

hoM%) = E n ^n(x), E n = 2n + 2k (8) 
it transforms them into the solutions of another Schrodinger equation 

<f n (x) = N n Lip n (x), hi(p n (x) = E n <f n (x) 

with the same eigenvalues E n . The factor N n is introduced to ensure that the states 
ifi n are normalized to unity. New exactly solvable Hamiltonian hi = ho + A(x) is 
defined by the potential difference A = A(x) = — 2(lnw)". The function u = u(x) 
called transformation functiorM is a solution to the initial Schrodinger equation 

h u(x) = au(x) 

with a ^ E . 

For the sake of simplicity we will restrict ourselves to the case when a < Eq, 
u(x) Va; G (0, oo) and l/u(x) is not a square integrable function in the interval 
[0, oo). In this case u ^ H and the set |J <f n )} constitutes a complete basis in the 
Hilbert space H provided the initial system {| ip n }} is complete.! In terms of the 
super symmetric quantum mechanics this case corresponds to a broken supersymme- 
try. 
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Let us now choose the following solution of the initial equation (||): 
u = u p {x) = x^-^Lf' 1 (y) exp (x 2 /a) , y = -x 2 /2, 

a = -2(k+p), p = 0,1,2,... 

as the transformation function. It is worth stressing that L 2k ^ 1 (y) 7^ Vx 7^ when 
y = —x 2 /2(< 0) and, hence, the function u p (x) is nodeless for x > 0. The function 
Lq(x) = L 0p (x) has in this case the form 

T ( \ - 1 ~ 4k X XL ?-1 ^ 
0p[X) 2x 2 Lf-^{y) 

and for the potential difference we obtain 

Ak- 1 



A = A p (x) = -l + 

x z 

Lf-, (y) 



x 2 L^(y) + L^. 1 (y) + ^ 



Lf- l (y) [Lf~ x (y)\ 



2 



-x 2 /2. 



It is worth noticing that because of the condition L 2 (y) 7^ for x 7^ the 

potential difference A p (x) has a single pole at x — 0. The asymptotic behavior of 

this function near the point x = is as follows: A p (x) — > —1 — p/k + (4k — l)/x 2 . 

In the limit as x — > 00 this function tends to —1. For 6 7^ the behavior of the 

transformed potentials V^(x) = Vo(x) + is nearly the same as the behavior 

of the initial potential Vq(x) = x 2 /A + b/x 2 . There is only a small difference in the 

region of the minimum. When 6 = the initial potential is regular at x = whereas 

all V p (x) are singular: V p (x) — > — 1 — Ap/3 + 2/x 2 . As an illustration, we give the 

explicit expression for two first potential differences 

a z \ , 4fc - 1 A , . . . . 4 32k 

A {x) = -1 + — , Ai(x)=Ao(x) + 



x 2 ' 1V y uv y 4A; + x 2 (4A; + x 2 ) i ' 

Thus, the transformed Hamiltonians hi = h^f 1 = —d 2 /dx 2 + V p (x) are bounded from 
below and are essentially self adjoint operators defined on the dense set in the Hilbert 
space H. All these Hamiltonians have only discrete spectrum. The position of the 
energy levels is the same as before transformation: E n = 2n + 2k, n — 0, 1, 2, . . .. 

The operator L + = —Lq(x) — d/dx realizes the transformation in the inverse 
direction 

|Vn> = N n L + \cp n ), (9) 
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and together with L participates in the following factorization: 

L + L = ho-a, LL + = h l -a . (10) 

We note that in our case operators L and L + are well-defined V0 G H and are 
conjugated to each other with respect to the inner product in the space H. This 
observation together with the factorization properties ( |10"D is very useful in calculat- 
ing integrals involving transformed functions. In particular, for the normalization 
constant N n we readily find N~ 2 = (ip n \hi — a\f n ) = E n — a = 2p + 4k + 2n. 

Now, by following the method of Ref.ffBI we find the transformed coherent states 



<p z {x) = Ni z Lip z (x) = b n z n <fn, 

n=0 

N = N 0z N lz /N , b n = a n N Q /N n 

From the factorization properties ( |T0"D one readily derives the value of the normal- 
ization constant 

, T _ 9 4A; + 2p — 2pzz 

N lz = 1 - • 

l-ZZ 

Using the identity operator 

oo 

Eb»)W = 1 ( n ) 

n=0 

one can calculate the integration measure v{z) which realizes the following identity 
decomposition: 

\<p g )(<p g \dv(z) = l. (12) 



\z\<l 

For this purpose we pass from the complex variables z and z to the polar ones z = 
y/xexpi(f) and look for the measure in the form dv = |/i (x) dxdcp. The integration 
over the angle variable <fr i n (P^2|) is trivial. Then making use of the formula ([11]), we 
obtain the equation for the function h{x) 



x n (1 - x) 2k+l n\Y (2k) 



7r / dxh(x] 



2k + p — px r (n + 2k) (n + 2k + p) 

So, we arrive to a moments problem well known in mathematics.il In our case 
the solution to this problem can readily be obtained with the help of the integral 
representation of the beta-function 

B ( a ,6)=£wif=/ I x-(i-,r^ 

1 [a + b) Jo 
9 
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and the following identity: 

r(n + l)r(2fc) * 2k -1 n , . , n . 

(n + 2k + p y(n + 2k) =^ 2k + P -j-l B ^ + ^^ 2k + p -^ 

where C J p are the binomial coefficients. It is not difficult to prove this identity by 
induction. The final expression for the function h(x) is 

2k - 1 p r-? M - r) p ~ j ~ 2 

Ml) = ^ (2t+p _^ ) g c ; |_(i_j_ . 

Once we know the measure we can construct the holomorphic representation 
(p(z) in the unit disk of a vector |</?) G H which is defined by its Fourier coefficients 
c n with respect to the base {|y? n )} 

oo 

(ipz\tp) = Nip (z) , if (z) = b n c n z n . 

n=0 

It follows from these relations the holomorphic representation for the transformed 
discrete basis vectors ip n (z) = b n z n and for the transformed coherent states 

With the help of the transformation operators L and L + we construct the lower- 
ing and raising operators p± = Lk±L + for the transformed discrete basis functions 

p + \ Vn ) = -(AT n AT ri+1 )- 1 [(n + \){n + 2k)} l ' 2 \ Vn+l ) , 

p_\ Vn ) = -(N n N n ^)- 1 [n(n + 2k - 1)] 1/2 ly^). 
Together with the operator p — they form a non-linear algebra 

\po,P±] = ±P±, \p-,p+] = 2 (2k(l - k) - p a + Apfj (2p - a) . 

It is interesting that the operators p± are of the fourth degree in derivative in the 
coordinate representation while in the holomorphic one they are of the second degree 

P„M = 4 + *, ,,.(*) = 2,^ + 2(2* +,,)£, 

d 2 ( d \ 

p+ (z) = 2z 3 — + 2z(2k + p + 2) iz— + 2k \ . 
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Since the operators L and L + realize the correspondence between the basis sets 
an d {Iv'n)} we can find their holomorphic representation 



d 



By assuming that the inner product of two vectors \(pi) and \<p2) is independent 
on the representation used and taking into account the identity decomposition (|T2|) 
we can equip the lineal £ with a new inner product 

((PiM = (Vi(z) \<Pz(z)) = / 7p 1 (z)(p 2 (z)e- f{1) diy(z). (13) 

J\z\<l 

It is naturally to declare the function 

/« = /«(*,*) = ln|fob.>r a = / (0) + ln (14 ) 

h — a 

a new Kahler potential. It then follows the expression for a new Hermit metric 

(i) = f (i) 2k + 1 p(2k + p) 

9 h ~ z {\-zzf {2k + p-pzzf 

The unit disk represents now a phase space with the varying Gauss curvature 
/C« = -2_g_i n ^(i). The explicit expression for K.^ is rather complicated so that 
be omitted here. We will note only that K,^ — > when k — > oo. 
The function 

(*,*) = f> n (*') = (1 - 2Jb + P - P ^ 

plays the role of the transformed Bergman kernel which takes off the integration 
with respect to the new inner product (|13|). 

By making use of the covariant Berezin symbols we can establish correspondence 
between the quantum mechanical operators and the classical observables which are 
the functions of a point in the unit disk. The function 

tt i it i \ a Akzz (2k + p + I - pzz) 

corresponds to the transformed Hamiltonian operator and therefore it is the trans- 
formed Hamilton function. Similarly, the functions P + and P_, 

P + = P_ = ( l p-\p + \ (p -) = 

11 
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Akz 



2k + p — pzz 



(2k + 1) (2k + 2) . iN 2fc+l 



^1 - zzY ' ' 1 — zz 

correspond to the lowering and raising operators. 

The change of the Kahler potential induces the change of the symplectic form 
which remains closed because of its dimension. As a consequence we have another 
Poisson brackets in the unit disk. The functions Pq = \H\ and P± form now neither 
a closed Lie algebra nor a polynomial one with respect to the new Poisson bracket 
but it is not difficult to check the following relations: 



{P ,P ± } = ±iP ± , {P ,z} = iz, {P ,z} 



-iz. 



Two latter relations mean that the equation of trajectories for the transformed clas- 
sical system in the transformed phase space remains unchanged under the Darboux 
transformation. 



VI. CONCLUSION 



Using the coherent states technique we have constructed a classical counterpart 
of the transformed a la Darboux quantum singular oscillator. It has been shown 
that the Darboux transformation for this system translates into the transformation 
of the Kahler potential. New Kahler potential has a supplementary summand of the 
form \n(H$ — a) where Hq is initial Hamilton function corresponding to the quantum 
Hamiltonian operator for the singular oscillator. This leads to the distortion of the 
initial phase space. However, the distortion is consistent with the change of the 
Hamilton function in such a manner that the curves of a constant energy (classical 
trajectories) in the transformed phase space remain unchanged. A similar obser- 
vation has recently been discusse d in Lagrangian formalism and for the canonical 
method of quantization.!!! 
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